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Abstract — We  present  a  studs'  of  scattering  from  objects  in  a  free 
unbounded  ocean  and  for  objects  near  an  interface  and  in  a  waveguide. 
Major  emphasis  is  on  results  and  discussion  of  theoretical  development  is 
limited  to  the  essentials.  An  examination  of  backscattered  echoes  is 
presented  for  elastic  targets.  Angular  distributions  are  examined  for  all 
target  types,  as  well  as  environments,  with  the  exception  of  objects  in  a 
waveguide.  In  a  waveguide,  beamforming  techniques  are  most  appropri¬ 
ate  for  detection  of  objects.  We  also  examine  rigid,  soft,  and  elastic 
objects  that  are  either  spheroidal  or  cylindrical.  Aspect  ratios  studied 
range  from  3:1  to  24:1,  and  kL/2  ranges  (*  is  wavenumber  and  L  is  the 
length  of  the  object)  are  from  3  to  120.  f  •'  | .  ■'  )  L  *  £  <  tur-.  ; 
s:  ■■  =>‘  7'/ 

'  I.  Introduction 

IN  classic  and  modem  physics,  the  scattering  of  waves  or 
particles  from  objects  or  interfaces  is  often  the  only  means 
by  which  their  properties  can  be  ascertained.  It  is  therefore  not 
surprising  that  scattering  theory  and  its  implementation  repre¬ 
sent  a  considerable  portion  of  the  articles  in  both  engineering 
and  physical  science  journals.  In  modem  physics  problems, 
one  deals  not  only  with  scattering  problems  but  also  with  the 
more  general  area  of  reaction  theory,  in  which  the  final 
products  of  some  physical  encounter  differ  from  those  of  the 
initial  quantities.  In  the  world  of  modem  physics,  relevant 
quantities  of  interest  may  include  angular  momentum,  nuclear 
radii,  and  nonclassic  quantities,  such  as  intrinsic  and  isospin, 
parity,  and  other  abstract  quantum  mechanical  entities.  Often, 
the  details  of  the  scattering  event  are  not  crucial  in  obtaining 
relevant  information,  and  simple  perturbation  theories  with 
vague  interactions  are  sufficient  to  yield  results.  Thus,  because 
of  the  broader  scope  of  abstract  quantities  that  enter  the 
problem,  modern  physics  scattering  and  related  reaction 
problems  are  more  difficult,  conceptually,  than  their  related 
classic  counterparts  but  they  are  often  simpler  from  the 
computational  point  of  view.  The  reason  for  this  difficulty  is 
that  classic  problems  involve  more  precise  knowledge,  such  as 
details  in  shape  and  boundary  conditions.  Classic  scattering  is, 
therefore,  more  concerned  with  shape  and  material  properties. 
Nevertheless,  a  common  set  of  concepts  connects  the  two 
areas  of  study  and,  ironically,  the  older  area  of  classic 
scattering  has  been  enriched  by  notions  and  techniques 
developed  in  the  more  recent  methodology  of  quantum 
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scattering  [1].  The  objective  of  this  research  is  to  develop 
classic  scattering  theory  in  a  consistent  manner  from  fairly 
simple  to  rather  complicated  physical  problems,  with  empha¬ 
sis  on  application  and  related  physical  interpretation.  The 
methodology  and  the  physical  ;nterpretation  of  the  work 
presented  here,  as  well  as  related  references,  draw  heavily 
from  modem  physics  concepts.  We  will  deal  exclusively  with 
objects  submerged  in  a  fluid  in  which  the  fluid  is  first 
unbounded,  then  near  an  interface,  and  finally  in  a  waveguide. 
The  objects  will  vary'  from  impenetrable  to  elastic  shells. 
Measurable  quantities  common  to  both  areas  of  physics  are 
monostatic  and  bistatic  angular  distributions  (differential 
scattering  cross  sections  in  quantum  physics)  and  backscat¬ 
tered  form  functions  (excitation  functions  in  quantum  physics) 
that  can  indicate  resonance  phenomena.  From  these  calcula¬ 
tions  we  will  show  how  it  is  possible  to  extract  shape  features 
and  material  properties  under  certain  conditions.  Emphasis  in 
this  study  is  on  application  and  interpretation  of  results,  with 
an  outline  of  the  theories  used.  Section  II  deals  with  the 
unbounded  problem.  Section  III  deals  with  scattering  from 
objects  near  an  interface.  Section  IV  deals  with  scattering  from 
objects  in  a  waveguide.  Section  V  summarizes  the  results. 

II.  Scattering  from  Submerged  Unbounded  Objects 

The  following  development  pertains  to  the  problem  of 
scattering  from  objects  submerged  in  a  fluid  in  which  there  are 
no  boundaries.  This  proves  to  be  a  good  approximation  for 
high-frequency  problems  in  deep-ocean  conditions  in  which 
the  object  is  far  from  any  boundary  and  detection  is  close 
enough  so  that  no  guided  waves  can  be  generated  from  the 
scattered  object.  In  subsequent  sections  we  will  consider 
problems  in  which  only  one  interface  is  significant  and, 
subsequently,  in  which  waveguides  that  are  either  refractive  or 
bottom  limited  are  important.  For  the  present  problem  one  can 
deal  with  a  variety  of  boundary  conditions  on  the  target.  The 
conditions  include  objects  that  can  be  impenetrable,  satisfying 
either  soft  (i.e.,  bubbles)  or  rigid  (very  heavy  objects) 
conditions.  Other  cases  can  involve  fluid  targets  and  elastic 
solids,  elastic  shells  in  which  the  shells  may  be  evacuated  or 
embedded  with  a  fluid  or  another  solid.  We  have  listed  these 
problems  in  order  of  difficulty  and  all  have  been  solved 
exactly  for  spheres  and  infinite  cylinders;  but  for  arbitrary 
shapes  the  problems  have  presented  theoretical  and  numerical 
problems.  We  indicate  how  these  problems  can  be  solved  by 
boundary  integral  methods  and  put  in  a  form  that  will  facilitate 
implementation  of  these  methods  in  the  more  complicated 
problems  involving  an  interface  or  a  waveguide.  The  form 
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used  draws  much  from  quantum  scattering  theory  and  was 
originally  developed  by  Waterman  [2]  for  electromagnetic 
problems  and  later  by  others  [3]-[7]  for  a  large  variety  of 
problems.  Improvements  to  the  method  have  also  been  made 
by  several  authors  and  recent  applications  of  considerable 
interest  have  been  presented.  The  following  outlines  the 
method. 

Let  us  suppose  that  we  have  an  incident  field  U‘  and  we 
wish  to  find  a  scattered  field  Vs.  Since  the  scattering  process  is 
linear,  we  can  obtain  the  scattered  field  from  the  incident  field 
by  means  of  a  linear  transformation,  which  can  be  obtained 
from  several  methods.  Since  this  operator  represents  a 
transition  from  the  incident  to  the  scattered  field,  it  is  usual  to 
denote  the  transformation  by  a  T  so  that  we  have 

(Jf=  TU‘.  (1) 

One  can  obtain  T  from  normal-mode  theory  for  spherical 
objects,  but  we  will  employ  the  extended  boundary  condition 
(EBC)  method  developed  by  Waterman  [3j,  since  it  works  for 
fairly  general  geometries  and  is  amenable  to  a  variety  of  useful 
algorithms.  The  EBC  method  is  often  called  the  T-matrix 
method,  but,  in  our  view,  misrepresents  the  method,  which  is 
really  a  boundary  integral  method. 

The  starting  point  for  the  EBC  method  for  fluids  is  the 
Helmholfz-Poincare  integral  representation  for  the  total  field 

U(r)  =  U'(r)  +  |  (U(r')dG(r,  r')/dn 

~G(r,  r')dU(r')/dn)  dS  (2) 

where  r  is  some  exterior  point.  If  r  is  some  interior  point,  then 
the  field  U  is  extinguished  so  one  has 

0=  t/'(r)+  j  ( U(r')dG(r ,  r')/dn 

-G(r,  r')dU(r')/dn)  dS.  (3) 

The  quantity  G  is  an  outgoing  Green’s  function,  and  r'  is 
associated  with  quantities  defined  on  the  surface  of  the  object. 
The  last  equation,  although  well  known,  is  usually  thought  to 
be  of  little  use.  However,  Waterman  [3]  recognized  that  it 
would  be  useful  as  a  constraint  to  eliminate  the  surface 
quantities  U(r')  in  (3)  by  solving  (3)  to  obtain  U(r')  to  be 
used  in  (2).  Note  that  we  have  two  equations  aVid  two 
unknowns,  namely,  the  scattered  field  and  the  surface  terms 
on  the  object.  Equation  (3)  can  be  solved  in  a  number  of  ways. 
It  can  be  solved  by  using  eigenvalue  techniques  [9],  [1 1]  or  by 
representing  it  in  terms  of  a  matrix  equation  by  expanding  the 
incident  field  and  the  Green’s  function  in  a  biorthogonal 
representation.  Both  the  Green's  function  and  the  incident 
field  are  given  below: 

f/'(r)  =  £  a,  Re  y>,(r)  (4a) 

G(r,  r')  =  iK  £  Re  *,-(#<)*,■(#■>).  (4b) 

The  quantities  Re  <p,  and  y?.  are  usually  incident  and  scattered 
partial  waves,  with  the  latter  satisfying  the  Sommerfeld 
radiation  condition.  We  still  do  not  know  the  surface  quantity 
U(r'),  but  we  can  expand  it  in  a  complete  set  defined  on  the 
surface  (i.e. ,  we  satisfy  closure).  Bv  doing  that  and  truncating 


at  some  point,  we  reduce  the  problem  to  that  of  a  matrix 
equation  in  which  we  can  identify  the  scattered  field  with  the 
incident  field.  We  will  not  include  the  details,  since  they  have 
been  presented  many  times.  We  will  give  only  detailed 
expressions  for  rigid  and  soft  cases,  but  not  for  the  more 
complicated  elastic  solid  targets  and  elastic  shells.  We  will, 
however,  present  calculations  for  the  elastic  case. 

Before  we  go  any  further  in  the  discussion,  we  must  find  a 
complete  set  of  functions  that  will  represent  the  surface 
quantities  U(r').  If  we  are  dealing  with  only  the  rigid  or  soft 
problem,  then  by  comparing  the  EBC  solution  with  the 
normal-mode  solution  we  can  show  that  they  are  formally 
identical  for  the  spherical  case  if  we  chose  the  following 
expansions: 

f/(r')  =  £  6,  Re  ¥>,(/-')  (5a) 

or 

VU(r’)  •  n  =  ^  bjV  Re  <pt(r')  •  A  (5b) 

where  (5a)  represents  the  expansion  for  the  rigid  case  and  (5b) 
represents  the  soft  case.  These  expressions  are  now  restricted 
to  spherical  representations.  Earlier  works  [3],  [8]  show  that 
these  functions  represent  a  complete  set  on  the  surface  for 
either  the  pressure  release  boundary  condition  or  the  rigid 
boundary  condition.  When  those  expansions  are  imposed,  one 
arrives  at  matrix  equations  that  relate  the  a’s  with  the  b’%  and 
the  /’ s,  where  the  /’ s  are  the  partial  wave  scattering 
coefficients  in  which  the  scattered  field  is  represented  by 

C /,<')«  £ /rtfc-  (6) 

The  y>/s  are  outgoing  partial  waves  as  stated  earlier  for  (4). 

For  the  case  of  a  three-dimensional  problem  in  a  fluid,  then 
the  <pi  s  are  just  outgoing  spherical  Hankel  functions  times 
spherical  harmonics.  For  the  incident  field,  the  quantity  Re  y>, 
is  the  spherical  Bessel  function  times  a  spherical  harmonic. 
Finally,  when  we  implement  (6),  in  (4)  and  (3)  we  will  get  two 
sets  of  equations,  which  are  represented  by  (7);  they  will  be 
the  matrix  equations  we  will  have  to  deal  with  to  find  the 
scattered  solution: 

/=  -  /  Re  Qb 

a=iQb.  (7) 

The  Q  matrices  appearing  previously  for  the  rigid  and  soft 
cases  are  given  in  (8): 

0,=  {Re»,^</S 

r  d  Re  <p. 

=  \'p,~diTdS-  (8) 

The  matrix  expressions  for  the  fluid  target  case  are  a  little 
more  complicated  but  are  not  much  more  difficult  to  calculate. 
One  can  see  that  (7)  can  be  used  to  eliminate  the  surface 
quantities  represented  by  the  b's.  Note  that  a,  b,  and  /  are 
column  vectors  that  are  coefficients  of  the  orthogonal  bases 
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functions  used.  We  can  eliminate  b  using  these  equations  to  first  examine  bistatic  angular  distributions  which  are  computed 
obtain  /  in  terms  of  a  to  obtain  from  the  following: 


/=  -Re  QQ  'a.  (9) 

From  the  definition  in  (1)  we  can  see  that  the  T-matrix  is 
given  by 


T=-ReQQ'.  (10) 

This  equation  must  satisfy  reciprocity  for  all  cases  and  the 
generalized  optical  theorem  for  nonlossy  environments.  These 
properties  can  be  checked  to  see  if  this  is  not  the  case.  This 
way  of  testing  determines  if  enough  terms  were  included  in  the 
partial  wave  expansions  represented  by  (4)— (6).  There  are 
other  ways  of  dealing  with  these  equations  in  which  one  can  be 
assured  of  obtaining  convergence,  and  there  are  methods  in 
which  this  has  been  dealt  with  as  discussed  by  Werby  and 
Chin-Bing  [9].  In  fact,  we  use  the  methods  in  that  paper,  in 
particular,  when  we  do  high-aspect-ratio  targets  that  will  be 
discussed  shortly.  In  the  following  results  the  coupled  high- 
order  T  matrix  had  to  be  used  for  the  problem  to  be  tractable. 
There  is  also  an  alternative  method  in  which  one  need  not  deal 
with  a  T matrix  as  such,  but /could  be  calculated  directly  from 
calculating  surface  quantities  by  an  eigenvalue  method  [10], 
[11].  We  wish  to  illustrate  calculations  for  solids  and  shells. 
We  will  only  indicate  the  forms  of  the  T  matrix  for  a  solid 
object  and  a  shell  in  (11)  and  (12),  where  (1 1)  is  given  in  [6] 
and  [12],  and  (12)  is  given  in  [13]: 

T=  -Re  QRP[QRP]  1  (11) 

T=  ~(QRR+QR0T2)M 'PKQor  +  QMM'P]-'  (12) 

where  M=RcR  +  RT2  +  iT2.  Here  T2  corresponds  to  the 
T  matrix  derived  from  scattering  from  the  inner  shell. 

The  definitions  of  the  matrices  for  R,  P,  and  the  various  Q’s 
are  defined  in  those  papers  and  represent  quantities  obtained 
by  the  additional  boundary  conditions  for  the  elastic  case, 
where  they  appear  as  constraints.  We  might  add  that  the  Q’s 
for  the  elastic  fluid  interface  cases  are  no  longer  square  but  are 
n  x  3n  matrices,  which  adds  some  complications  to  the 
problem  that  are  somewhat  more  difficult  to  deal  with.  Werby 
et  at.  [1 1],  [14]  indicate  how  some  of  these  equations  can  be 
dealt  with  in  a  more  effective  manner.  A  transformation 
method  was  also  developed  [11]  to  show  how  some  of  the 
inherent  numerical  difficulties  obtained  by  these  expressions 
can  be  overcome.  Our  aim  here  is  to  present  results  that  deal 
with  fairly  complicated  problems.  In  the  next  section,  mainly 
numerical  results  of  scattering  from  unbounded  objects  will  be 
presented  for  very  high-aspect-ratio  rigid  and  soft  problems. 
Elastic  solids  and  shells  for  some  of  the  interesting  effects  due 
to  resonances  will  also  be  examined. 

For  the  first  example  we  chose  a  150-m-long  spheroid  that 
has  a  semiminor  axis  of  5  m  (aspect  ratio  of  15:1).  Here,  we 
chose  kL/2  =  120,  where  k  is  the  wavenumber  and  L  is  the 
length  of  the  object.  We  examine  bistatic  angular  distributions 
for  several  cases.  The  incident  field  is  chosen  to  be  a  plane 
wave  that  will  impinge  on  the  object  at  the  following  angles: 
0°,  30°,  60°,  and  90°  relative  to  the  axis  of  symmetry.  We 


m= 


2  AnPn  (C0S  ( 9))T" 


(13) 


where  6  is  the  angle  relative  to  the  axis  of  symmetry.  Note  also 
that  the  A 's  are  dependent  on  the  azimuthal  index  M  and  are  a 
function  of  the  orientation  of  the  incident  field  relative  to  the 
axis  of  symmetry  of  the  object. 

Fig.  1(a)  illustrates  the  case  for  which  the  incident  field 
enters  the  object  along  the  axis  of  symmetry.  At  this  frequency 
the  field  is  altered  predominantly  in  the  forward  direction.  Fig. 
1(b)  and  (c)  is  more  interesting,  since  it  illustrates  a  principle 
similar  10  Snell’s  law  for  a  flat  surface.  Here,  we  have  a 
diffracted  wave  in  the  direction  of  the  incident  field  and  also  a 
reflected  wave.  Because  the  object  has  such  a  high  aspect 
ratio,  it  behaves  almost  as  a  flat  surface  for  a  large  part  of  the 
surface  area;  thus  the  angle  of  incidence  and  the  angle  of 
reflection  are  the  same.  Fig.  1(d)  illustrates  the  case  for 
scattering  perpendicular  to  the  axis  of  symmetry,  or  broadside. 
Note  that  at  this  kL/2  the  scattered  field  is  well  focused  and 
has  the  largest  amplitude  in  the  forward  direction.  It  has  a 
target  strength  that  is  26  dB  above  the  backscattered  contribu¬ 
tion.  This  difference  is  independent  of  the  actual  size  of  the 
object  and  is  only  a  function  of  the  shape  and  the  quantity  kL/ 
2.  The  calculations  in  Fig.  I  (a)— (d)  were  done  for  a  fixed 
frequency.  By  the  time  kL/2  =  120  the  scattered  field  is  on 
the  verge  of  approaching  the  geometric  limit.  For  lower 
frequencies  the  pattern  would  be  more  wave-like.  To  see  this 
effect  let  us  examine  end-on  and  broadside  scattering  for  a  75- 
m-long  spheroid  with  a  semiminor  axis  of  5  m.  Fig.  2 
illustrates  the  case  for  kL/2  =  7.5,  15,  30,  60,  and  120.  The 
top  curves  are  end-on  and  the  bottom  curves  are  broadside. 
Note  that  the  figure  illustrates  wave-like  properties  at  the  lower 
frequencies  but  becomes  more  directional  and  predominate  in 
the  forward  direction  with  increasing  frequency.  Moreover, 
the  forward  scattering  increases  with  increasing  kL/2,  since 
most  of  the  field  goes  off  in  the  forward  direction  for  the 
higher  kL/2's.  Fig.  3  illustrates  this  point,  as  we  plot  the 
forward  target  strength  for  broadside  scattering  as  a  function 
of  frequency  that  corresponds  to  the  kL/ 2's  in  Fig.  2.  Note 
that  the  quantity  grows  rapidly  at  first  but  appears  to  level  off 
at  the  higher  frequency  end.  Our  plots  are  presented  as  a 
function  of  kL/2  in  which  we  add  the  geometrical  cross 
section  to  the  form  function.  If  we  were  to  plot  the  form 
functions  as  a  function  of  ka  for  the  broadside  case,  where  a  is 
the  semiminor  axis,  and  if  we  were  to  exclude  the  geometrical 
cross  section,  then  for  objects  of  a  fairly  high  aspect  ratio,  the 
above  quantity  would  be  about  the  same  in  magnitude.  That 
means  the  forward-scattering  form  function  of  an  object  with 
an  aspect  ratio  of  7.5  and  kL/2  =  60  would  have  about  the 
same  magnitude  for  broadside  scattering  as  one  with  an  aspect 
ratio  of  15  but  kL/2  =  120,  since  both  have  ka  =  8.  This 
implies  that  high-aspect-ratio  objects,  when  scattered  from 
broadside,  behave  somewhat  like  infinite  cylinders,  as  we 
might  expect. 

Turning  to  monostatic  angular  distributions,  we  find  these 
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quantities  to  correspond  to  measuring  the  backscattered  form 
functions  (echos)  as  the  receiver  rotates  360°  about  the  object. 
For  the  high-frequency  case,  this  would  give  an  outline  of  the 
object  (its  shadow),  since  we  would  be  close  to  the  optical 
limit.  However,  for  the  low-frequency  case,  this  would  yield  a 
diffraction  effect  due  to  interference  between  the  waves 
reflected  from  different  points  of  the  object.  We  have  shown 
that  the  location  of  the  nulls  (as  exhibited  in  Fig.  4)  and  the 
number  of  nulls  are  almost  an  exclusive  function  of  the  kL/2 
value.  This  means  that  a  knowledge  of  the  number  or  location 
of  the  nulls  for  elongated  objects  can  determine  the  length  of 
the  object,  since  the  equations  below  relate  kL/2  in  a  1:1 
correspondence  between  the  number  of  nulls  and  the  location 
of  the  nulls  [20J: 


6„  =  arcsin 
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Fig.  I.  Scattering  from  a  rigid  spheroid  of  aspect  ratio  15:1  for  kL/2  = 
120:  (a)  end-on  incidence:  (b)  30°  relative  to  the  axis  of  symmetry;  (c)  60° 
relative  to  the  axis  of  symmetry;  and  (d)  90°  relative  to  the  axis  of 
symmetry. 
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The  top  three  plots  in  Fig.  4  were  performed  for  objects 
with  aspect  ratios  of  24: 1  and  for  kL/2  =  40.  The  First  of  the 
top  three  objects  is  a  rigid  spheroid,  the  second  is  a  soft 
spheroid,  and  the  third  is  a  rigid  cylinder.  Note  that  the 
patterns  are  quite  different  for  the  three  cases.  The  soft  object 
exhibits  little  diffraction,  which  is  a  well-known  fact  from 
backscattered  form  functions  for  soft  objects  [15),  [16[.  In 
fact,  plots  of  bistatic  angular  distribution  of  soft  objects  are 
much  more  focused  than  their  rigid  counterparts  due  to  less 
diffraction  effects,  which  is  illustrated  in  Fig.  5.  The  second 
and  third  sets  of  objects  correspond  to  monostatic  angular 
distributions  for  rigid  spheroids  for  aspect  ratios  of  15:1  and 
7.5:1  for  kL/2  values  from  left  to  right  of  30,  60,  and  120, 
respectively.  Note  that  the  patterns  approach  the  optical  limit 
with  increasing  frequency.  The  rate  of  approach  is  more  rapid 
for  the  lower  aspect  ratio  case,  since  the  ka  values  for  this  case 
are  double  that  of  the  higher  aspect  ratio  case. 

Fig.  5  illustrates  the  case  of  soft  scatterers  for  bistatic 
angular  distributions  for  end-on,  30°,  60°,  and  90°  relative  to 
the  axis  of  symmetry  for  an  aspect  ratio  of  15:1  for  kL/2  = 
30.  Note  that  due  to  decreased  diffraction  effects  for  soft 
objects,  this  scatterer  behaves  less  wave-like  than  its  rigid 
counterpart  and  tends  to  focus  the  field  more  sharply. 

Backscattered  form  functions  are  target  echoes  plotted  as  a 
function  of  frequency.  Only  elastic  and  fluid  targets  sub¬ 
merged  in  water  can  support  resonances.  (We  exclude  the  case 
of  bubbles.)  We  illustrate  this  case  by  examining  end-on 
incidence  of  a  steel  elastic  spheroid  where  the  compression;!! 
velocity  is  5.95  km/s,  the  shear  velocity  is  3.24  km/s,  and  the 
density  is  7.7  relative  to  that  of  water.  The  plot  in  Fig.  6  is 
from  kL/2  =  4-16.  We  observ  two  interesting  features.  The 
First  feature  is  the  periodic  behavior  known  as  Franz  wav  es 
They  correspond  to  an  interference  between  the  specular 
waves  and  the  circumferentially  diffracted  waves.  The  mini¬ 
mum  and  maximum  values  correspond  to  destructive  and 
constructive  interference.  This  effect  is  strictly  geometrical 
and  would  be  identical  to  that  of  a  rigid  object  of  the  same 
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Fig.  2.  Scattering  along  the  axis  of  symmetry  and  broadside  for  a  spheroid 
75  m  long  and  10  m  wide  for  kL/2  =  7.5,  15,  30,  60,  and  120. 
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Fig.  5.  Scattering  from  a  soft  spheroid  150  x  10  m  for  kL/2  =  30:  (a) 
along  the  axis  of  symmetry;  (b)  30”  relative  to  the  axis  of  symmetry;  (c)  60” 
relative  to  the  axis  of  symmetry;  and  (d)  90”  relative  to  the  axis  of 
symmetry . 


shape.  The  difference  between  rigid  and  elastic  solids  arises 
due  to  the  coupling  of  transverse  and  longitudinal  waves  at  the 
fluid-solid  interface  [17).  These  result  in  interface  waves  and, 
for  this  case,  they  arc  leaky-type  Rayleigh  waves  that  result  in 
resonances  at  discrete  frequencies.  This  resonance  effect  is 
characterized  by  a  minimum  followed  by  a  sharp  rise,  which 
reflects  the  fact  that  these  types  of  resonances  change  phase 
rapidly  by  180°  relative  to  the  nonresonance  background  and. 
thus,  are  out  of  phase  and  then  rapidly  in  phase  with  the 
background  over  a  short  change  in  kL/2. 

Lower  aspect  ratio  targets  and  spheres  usually  manifest  a 
much  more  pronounced  null,  but  our  studies  indicate  that  with 
increasing  aspect  ratio  this  effect  is  less  evident  for  end-on 
resonances  due  to  the  already  low  value  of  the  specular 
reflection.  Note  that  there  are  two  Rayleigh  resonances  here, 
the  second  being  weaker  than  the  first. 

For  completeness,  we  illustrate  the  resonance  response  for  a 
very  thin  steel  shell  with  an  aspect  ratio  of  3:1  in  Fig.  7.  The 
parameters  for  steel  are  the  same  as  for  the  case  above.  We 
plot  the  backscattered  form  function  for  kL/2  from  3  to  15. 
Here  we  observe  resonances  at  two  locationc  These  reso¬ 


nances  correspond  to  modal  vibrations  that  result  in  standing 
waves  at  discrete  frequencies:  a  simple  picture  of  this  behavior 
has  been  demonstrated  [18].  Note  that  these  resonances  are 
also  characterized  by  a  minimum  followed  by  a  sharp  rise.  It 
has  been  demonstrated  that  the  proper  background  for  very 
thin  shells  is  sound  soft,  and  that  is  the  case  for  this  calculation 
[15],  [16],  [19], 

We  now  investigate  scattering  from  objects  near  an  inter¬ 
face.  Since  the  theory  on  this  subject  has  not  been  well  covered 
in  the  literature,  we  present  a  more  extensive  theoretical 
development  of  this  topic  than  for  the  other  cases. 

III.  Acoustic  Scattering  from  an  Object  in  a  Half-Space: 

The  Exact  Image  Solution 

The  calculation  of  the  scattered  acoustic  field  from  a 
nonspherical  object  in  free  space  is  not  a  simple  task.  This 
problem,  however,  has  been  successfully  treated  using  the  T- 
matrix  method.  The  T-matrix  method  is  due  to  Waterman  [2]. 
[3]  and  has  been  developed  for  high-aspect-ratio  objects  by 
Werby  and  Chin-Bing  [9]  and  Werby  [20].  The  scattering 
problem  is  even  more  difficult  when  the  object  is  near  a  plane 
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Fig.  7.  Form  function  of  steel  spheroidal  shell  of  aspect  ratio  3: 1  from 
AX/2  =  3-15. 

boundary.  Our  approach  is  to  use  the  free -space  T-matrix 
scattering  theory  and  apply  it  to  an  object  in  a  half-space  via 
the  method  of  images. 

The  free-space  scattering  theory  to  be  used  is  from  Peterson 
and  Strom  (21]  and  involves  the  T-matrix  formulation  for 
acoustic  scattering  from  an  arbitrary  number  of  scatterers.  The 
object  near  a  plane  boundary  is  replaced  by  the  object  in  free 
space,  along  with  its  image  and  an  appropriate  restriction  on 
the  incident  wave.  We  then  apply  the  theory  of  Peterson  and 
Strom  [21]  to  the  combined  system  of  two  objects.  The 


scattered  field,  restricted  to  the  half-space,  is  the  same  as  the 
scattered  field  from  the  single  object  in  the  half-space. 

Consider  now  an  object  in  a  half-space  bounded  by  a  plane 
surface  (z  -  0)  as  depicted  by  Fig.  8.  The  object  is  assumed  to 
be  rigid  for  the  sake  of  definiteness.  Only  the  free-space  T 
matrix  of  the  object  is  needed  so  other  types  of  objects  may  be 
used  in  the  same  formulation.  The  plane  surface  is  assumed  to 
be  rigid.  The  case  when  the  plane  surface  is  free  will  be 
discussed  later.  The  surface  of  the  object  is  denoted  by  St.  The 
point  0,  is  in  the  interior  of  the  object  at  a  distance  d  above  the 
plane  surface.  We  assume  that  the  object  is  insonified  by  a 
known  incident  time-harmonic  acoustic  wave  in  the  half¬ 
space.  The  total  field  is  the  sum  of  the  incident  field  and  an 
unknown  scattered  field: 

M(r)  =  «,(r)  +  «J(r).  (15) 
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The  vector  r  =  (x,  y,  z)  is  measured  from  the  point  0  on  the 
plane  surface  below  0|.  The  total  field  satisfies  the  Helmholtz 
formula 

U((p  =  u,(r)+  jj  u  +  (7')VgMf(7,  ?')  •  f  ds'  (16) 

for  7 in  the  half-space  outside  Si  and  7 inside  Si,  where  u .  is 
the  limit  from  the  outside  of  the  total  field  u  and  rj  is  the 
normal  to  the  surface  S|.  The  Green’s  function  for  a  half-space 
is 

gt*Ar,  r')=g(r,  7')+g(7,  R") 

where 

7'  =  (*',  y' ,  z') 

7"=(x',y\  - z ') 
and 

g(r,  7')  =  exp  (ik\  r-7'\)/(4x\r-7'\) 

is  the  free-space  Green’s  function. 

The  integral  equation  in  (16)  will  now  be  extended  to 
include  the  region  below  the  plane  z  =  0.  Assume  that  the 
incident  field  and  the  Green’s  function  for  the  half-space  are 
extended  as  even  functions  below  the  plane  z  =  0.  Let  S2  be 
the  image  of  S,  below  the  plane  z  =  0  as  shown  in  the  cross 
section  in  Fig.  9.  The  total  field  is  extended  as  an  even 
function  to  the  region  outside  S2  below  the  plane  2  =  0.  The 
right  side  of  ( 16)  defines  a  function  that  is  even  about  the  plane 
2  =  0.  This  function  is  the  same  as  the  total  field  in  the  region 
outside  S2  below  the  plane  2  =  0  and  must  be  zero  inside  S2. 
Therefore,  the  total  field  satisfies  the  integral  equation 

U^p  =  Uj(7)+  jj  u  +  (7')Vg(r,  r')  ■  rj  ds' 

’i 

+  jj  ut(rmWg(r,7")  vds”  (17) 

'2 

for  7 outside  S\  and  S2  and  7 inside  St  or  S2  Th?  integral  over 
S|  has  been  split  into  two  terms  containing  the  free-space 
Green’s  function,  and  the  variable  of  integration  in  the  second 
term  has  been  changed  to  7"  using  u,(r")  =  ut(r')  The 
integral  equation  in  (17)  is  the  same  as  the  one  that  would  be 
obtained  using  the  method  of  Peterson  and  Strom  [2 1  ]  for  two 
identical  rigid  objects.  If  the  plane  surface  is  a  free  surface, 
then  the  incident  wave  and  the  half-space  Green’s  function  are 
extended  as  odd  functions.  The  same  integral  equation  is 
obtained  except  that  the  total  field  is  odd  about  the  plane  2  =  0 
in  the  region  outside  5|  and  S2. 

Based  on  this  example,  we  proceed  by  using  the  method  of 
images.  We  follow  the  development  of  Peterson  and  Strom 
[21]  for  two  scatterers  to  solve  the  problem  of  scattering  from 
a  rigid  object  in  a  half-space.  The  incident  field  will  be 
expanded  in  regular  spherical  waves  relative  to  the  origin  0  as 


Fig.  9.  Schematic  of  objeci  and  its  image 

The  scattered  field,  far  away  from  0,  is  assumed  to  be 
expanded  in  spherical  waves  as  follows: 

Us(7)  =  £  fn<Pn(7). 

n 

The  goal  is  to  find  the  T  matrix  7X1,  2).  which  relates  the 
unknown  f„  to  the  known  a„. 

The  total  field  satisfies  (17),  which  can  be  written  more 
compactly  as 

U(j^  =  u,(7)  +  jj  uAr')Vg(ryr')  ■  n  ds'  (18) 

V*  2 

for  7 outside  S,  and  S2  and  7  inside  S,  or  S2,  where  T  varies 
over  the  surface  of  the  object  S|  and  the  surface  of  the  image 
S2.  All  the  functions  appearing  in  (18)  will  be  expanded  in 
spherical  waves.  This  expansion  yields  three  infinite  systems 
of  linear  equations  resulting  from  7  inside  S|,  rinside  S2,  and  7 
outside  S|  and  S2.  In  matrix-vector  notation  these  systems  of 
equations  are 

r(  -  3)a  =  iQ5\  +  ia(  -  23)  Re  Qa2  (19) 

R(ct)a  =  /a(2(7)  Re  £>a,  +  iQa2  (20) 

/=  -  iR(3)  Re  Qa ,  - 1 R(  -  3)  Re  Qa2.  (21) 

The  vector  3  =  (0.  0.  d)  is  the  position  vector  of  0,  and  -  3 
=  (0,  0,  -  d)  is  the  position  vector  of  02.  The  vectors  a,  and 
a2  contain  the  expansion  coefficients  of  the  surface  fields  on  S, 
and  S2.  The  matrix  Q  is  defined  by 

Qn  m=  -k  |  j  Re  ipn(r[)Vspn(r[)  •  rf  ds ' 

vi 

where  r'  is  a  position  vector  of  a  point  on  the  surface  S| 
measured  from  0,.  The  matrix  Re  Q  is  obtained  by  taking 
regular  functions  in  both  factors  in  the  integrand.  For  two 
scatterers,  two  matrices  Q ,  and  Q2  would  normally  be  needed; 
but  since  both  scatterers  are  identical,  the  common  value  is 
called  Q.  The  matrices  R  and  o  are  determined  by  the 
translation  properties  of  Re  \pn  and  and  can  be  found  in 
Peterson  and  Strom  [21]. 

The  T  matrix  7X1,  2)  which  satisfies 

7=7X1,  2)3 


“/(»')  =  2  a"  Re  *»( 7 )■ 


WERBY  AND  EVANS:  SCATTERING  FROM  OBJECTS  SUBMERGED  IN  OCEANS 

is  found  by  eliminating  and  a-,  from  (19)-(21).  If  bx  = 
iQa  i  and  b2  =  iQa 'hen  (19)  and  (20)  become 

. . v  - 3)a =  5|  -  <j(  - 2 3)Tb2 
R(3)a  =  -o(23)TBt  +  b2 

where  T  =  -  Re  QQ  1  is  the  free-space  T  matrix  for  the 
object  in  a  spherical  wave  representation.  Solving  for  h,  abd  b2 
yields 

£,  =  (1  -o(-23)T(2d)T)  x[R(-3)  +  a(-23)TR(3))a 

(22) 

b2  =  (1  -  a(2 3)To(  -2 d)T)  '\R(3)  +  o(23)TR(  -  3)\a. 

(23) 

The  expansion  coefficients  for  the  scattered  field  are  obtained 
from 

J=R(d)Tbt+R(-d)Tb2 ■  (24) 

Equations  (22)-(24)  reduce  to  the  formula  for  7~(1 ,  2)  in  (21, 
eq.  (33)| . 

For  the  specific  problems  to  be  considered,  the  incident 
field  will  be  taken  as  a  sum  of  plane  waves  that  satisfy  the 
boundary  condition  on  the  plane  surface  as  follows: 

u,(r~)  =  exp  (ik  ■  r)  +  exp  (ik  ’  ■  T) 

where  k  =  (k  cos  9,  0,  k  sin  9)  and  k‘  =  (k  cos  6,0,  -  k  sin 
0).  The  +  sign  is  used  for  the  rigid  plane  surface  and  the  - 
sign  is  used  for  the  free-plane  surface. 

We  now  examine  two  examples  in  which  we  scatter  from 
objects  near  a  rigid  interface.  Let  us  choose  the  case  of  a 
spheroid  some  distance  from  a  rigid  bottom  in  which  the  axis 
of  symmetry  of  the  spheroid  is  perpendicular  to  the  interface. 
Fig.  10  is  a  schematic  of  the  problem  under  investigation.  The 
solid  arrows  indicate  the  features  seen  in  the  absence  of  a 
bottom,  and  the  dashed  arrows  show  the  features  that  appear 
because  of  the  presence  of  the  bottom.  Please  note  that  all 
measurements  are  made  in  the  physical  region  and  do  not 
include  the  nonphysical  region  across  the  rigid  boundary.  If 
the  incident  field  is  at  an  angle  (relative  to  the  bottom)  greater 
than  zero  but  less  than  90°.  then  a  four-lobed  figure  would  be 
expected.  One  lobe  corresponds  to  a  reflection  from  the 
object;  a  second  lobe  corresponds  to  the  field  diffracted  from 
the  object.  Of  the  last  two  lobes  (which  arise  only  in  the 
presence  of  the  bottom),  one  corresponds  to  a  reflection  of  the 
reflected  wave  from  the  object  (which  would  come  off  in  the 
opposite  direction  of  the  incident  field),  and  the  other 
corresponds  to  a  reflection  of  the  diffracted  wave  from  the 
bottom. 

We  will  examine  these  cases  for  two  examples.  The  distance 
from  the  surface  is  eight  times  the  diameter  of  the  object.  In 
the  first  example  we  will  give  comparisons  of  the  object  near 
the  interface  with  one  in  a  free  environment.  The  object  is  a 
rigid  spheroid  of  aspect  ratio  2  and  kL/2  -  10.  Fig.  1 1(a) 
illustrates  an  example  in  w  hich  we  scatter  perpendicular  to  the 
bottom  (end-on  for  the  object).  For  the  free-space  example. 
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Fig.  10  Expected  response  from  scattering  of  object  near  a  half-space. 

little  backscattering  occurs,  and  most  of  the  scattering  is  in  the 
forward  direction.  However,  for  the  half-space  problem,  the 
results  are  characterized  by  a  reflection  from  the  rigid  bottom. 
These  results  look  much  like  what  is  obtained  when  scattering  is 
perpendicular  to  a  rigid  bottom.  Fig.  1 1(b)  and  (c)  illustrates 
the  examples  for  scattering  at  angles  of  30°  and  60°, 
respectively.  Note  that,  as  expected,  we  obtain  four  lobes  for 
the  half-space  problem  and  only  two  for  the  free-space 
problem.  Finally,  when  scattering  perpendicular  to  the  axis  of 
symmetry  of  the  object  and  parallel  with  the  bottom,  we  see  an 
almost  symmetric  distribution  for  the  half-space  problem, 
whereas  for  the  free-space  problem  most  of  the  scattering  goes 
off  in  the  forward  direction. 

It  is  interesting  to  perform  this  calculation  for  a  more 
elongated  object  and  at  a  higher  frequency,  since  the  effect  of 
the  reflected  waves  would  be  enhanced.  Fig.  12  is  an  example 
for  a  spheroid  of  aspect  ratio  10. 1  for  kL/2  =  30.  Here  Fig. 
12(a)-(d)  corresponds  to  scattering  from  0°,  30°,  60°,  and 
90°  relative  to  the  axis  of  symmetry  for  only  the  half-space 
problem.  The  reflected  and  diffracted  components  are  more 
pronounced  and  better  defined,  as  expected  for  the  higher 
frequency  cases.  From  these  examples,  objects  in  the  presence 
of  a  boundary  clearly  will  have  angular  distributions  quite 
different  from  those  in  a  free  fluid  or  far  from  a  boundary. 

IV.  Scattering  in  a  Waveguide 

Scattering  from  an  object  in  a  waveguide  becomes  a  very 
difficult  problem  if  one  attempts  to  treat  it  exactly.  Not  only  is 
it  difficult  to  formulate  for  equations  suitable  for  computation, 
but  even  approximate  methods  resulting  from  an  exact 
framework  result  in  expressions  that  are  not  very  useful  from 
the  caiculational  point  of  view.  We  therefore  resort  to 
approximations  based  on  physical  insight. 

To  do  this  let  us  take  an  object  in  a  waveguide.  Let  us 
assume  that  the  object  is  insonified  by  the  guided  wave  and 
that  it  is  small  enough  that  it  will  not  affect  the  incoming  field 
significantly.  Then,  to  obtain  the  scattered  field  in  a  region 
near  the  object  but  far  enough  away  from  the  boundaries  so 
that  it  will  not  have  a  chance  to  form  a  guided  wave,  one  can 
construct  the  scattered  field  using  the  T  matrix  in  conjunction 
with  the  incident  waveguide  by  representing  the  incident  field 
in  a  spherical  representation  so  as  to  be  compatible  with  the  T 
matrix,  since  the  T  matrix  is  in  that  representation.  One  can 
use  the  scattered  field  in  the  local  vicinity  of  the  object  to  be 
the  source  term  of  a  Hemholtz  equation  that  will  form  a  new 
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guided  wave  as  the  field  begins  to  interact  with  the  boundary . 
The  inhomogeneous  equation  is 

V2U+K2U=f  (25) 

where 

/=  Ta. 

Here,  we  construct  a  from  normal-mode  solutions,  which 
can  be  expressed  in  a  representation  of  incoming  and  outgoing 
plane  waves.  The  most  general  case  would  require  that  the 
individual  components  (pairs  of  incoming  and  outgoing  plane 
waves)  enter  the  target  with  different  wavenumbers  and  at 
different  angles  relative  to  the  target.  This  would  entail 
constructing  a  T  matrix  for  each  wavenumber.  One  would  then 
represent  the  individual  plane-wave  components  in  the  form  of 
a  Rayleigh  series  that  consists  of  products  of  spherical  Bessel 
functions  and  spherical  harmonics  that  yield  the  desired 
spherical  representation  of  the  sum  of  the  plane  wave.  One  is 
then  in  a  position  to  operate  on  each  T  matrix  to  obtain  the 
scattered  field  components.  These  components  are  then 
appropriately  summed  to  obtain  the  final  field,  which  is 
accurate  in  the  vicinity  of  the  field.  Once  /  is  known,  the 
solution  to  this  problem  can  be  obtained  in  one  of  several 
ways.  The  volume  integral  representation  is  as  follows: 

U'(r)=  \  G(r,  r')f(r')dr'  (26) 

where  G  is  the  approximate  Green's  function  obtained  from 
the  waveguide.  The  reason  that  we  do  not  integrate  over  the 
surfaces  of  the  upper  and  lower  bounds  as  required  in  other 
methods  is  that  the  Green's  function  here  is  associated  with  the 
waveguide  and  is  not  the  outgoing  Green's  function.  All  the 
information  contained  at  the  boundaries  is  thus  built  into  the 
waveguide  Green’s  function.  The  above  volume  integral  is  not 
very  efficient  for  calculation,  so  we  will  construct  two 
alternatives.  The  first  involves  a  surface  integral  representa¬ 


tion  of  the  problem,  namely 

t/s(0=  j  (f(r’)dG(r,  r')/8n-G(r ,  r’)df(r')/dn)  dS 

(27) 

where  the  surface  is  chosen  at  a  suitable  region  circumscribing 
the  object.  Alternative  techniques  to  the  above  method  are  one 
in  which  the  field  scattered  from  the  object  is  matched  at  some 
point  to  a  form  of  the  waveguide  solution  via  logarithmic 
derivatives,  or  one  in  which  certain  orthogonal  properties  of 
normal-mode  series  are  exploited  to  project  the  scattered 
solution  onto  the  normal-mode  representation.  Some  weighed 
residual  method  in  which  the  test  functions  are  normal-mode 
functions  can  also  be  used.  Currently  we  are  examining  these 
methods  for  comparison,  and  the  results  will  be  reported  in  a 
future  work.  The  following  is  a  development  of  the  method  we 
deem  most  suitable  to  solve  this  problem.  As  discussed  above, 
we  determine  the  field /just  outside  the  object.  In  particular, 
the  field  about  a  sphere  circumscribing  the  object  (suitably 
near  the  object)  with  origin  at  the  center  of  object  can  be 
obtained  from  the  expression 

/(p,  e,  *)=s  2  ^n,Tz,MY:(0Kn,  <t>nn) 

nm  n 

•  Y”,{dr,,  4,r.)h(Knp)/{K„pY'2  (28) 

where  the  a’”n , ,  are  projection  coefficients  of  the  normal¬ 
mode  functions  onto  the  spherical  (partial  wave)  solutions. 

For  the  unperturbed  waveguide  (70  in  normal-mode  repre¬ 
sentation  we  have 

U0=  \/2(e'*  *)/d  2  6n(y»zs)4>„{ynz)e‘«r/\jKHr.  (29) 

n 

The  Green's  function  is 

G0  =  (i/2d)  2  <t>n(ynZs)<t>n(ynZ)H0{K„r)  (30) 
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Fig.  12.  Scattering  from  an  object  of  aspect  ratio  10:1  for  kL/2  =  30:  (a) 
end-on;  (b)  30°  relative  to  the  plane;  (c)  60°  relative  to  the  plane;  and  (d) 
90°  relative  to  the  plane. 


when  H0  is  an  outgoing  cylindrical  Hankel  function  and  <t>„,  y,  where 

and  k  are  the  normal-mode  functions  and  the  vertical  and  f  ,vl/, 

horizontal  eigenfunctions,  respectively.  a„- y[f(p  ,  0,  <t>)b<t>„(y„z  )e*n  /(K„r  )  /bn 

If  we  insert  (28)  and  (30)  into  (27)  this  leads  to  an 

expression  far  from  the  object  of  the  form  ~  <t>n(y*z  )e’"r  V(xnr  (0))  bf(p  ,  6,  <P)/bn]  dS.  (32) 


Us(z,  r)  =  2  anWMynZte'n'/fary'2 


(31)  Note  that  (31)  is  in  the  form  of  a  normal-mode  solution,  and 
therefore  the  scattered  wave  forms  a  guided  wave  suitably  far 
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Fig.  13.  Beam-formed  response  of  rigid  spheroid  100  m  long  and  10  m  wide 
as  i(  passes  midway  between  a  source  and  receiver. 


from  the  object,  as  one  would  expect.  This  solution,  in  fact,  is 
continuous  throughout  space  and  satisfies  the  boundary  condi¬ 
tions  of  the  confined  environment  that  forms  the  guided  wave. 
We  will  illustrate  an  example  in  which  we  use  an  approximate 
solution  due  to  the  Huyghen  methods  described  above. 

We  will  now  present  a  calculation  in  which  we  scatter  from 
a  rigid  object  with  an  aspect  ratio  of  10: 1 .  The  object  chosen 
for  study  is  a  spheroid  that  is  100  x  10  m  in  a  water  column  of 
150  m.  A  simple  fluid  bottom  with  a  smooth  ocean  surface  is 
assumed.  The  speed  of  sound  in  the  water  is  chosen  to  be  1500 
m/s  and  the  frequency  is  250  Hz.  All  the  former  methods 
allow  for  a  variable  velocity  profile,  but  for  simplicity  we 
choose  an  isovelocity  profile  for  this  example. 

For  the  waveguide  the  bottom  is  chosen  to  be  a  sand  half¬ 
space  with  a  compressional  velocity  of  1733  m/s,  a  density  of 
1.95  relative  to  water,  and  an  attenuation  constant  of  0.97  dB 
per  wavelength.  The  source  and  receiver  are  situated  10  km 
apart,  and  we  allow  the  object  to  pass  midway  between  the 
source  and  the  receiver.  The  source  has  an  intensity  of  166  dB 
and  the  receiver  is  a  16-hydrophone  vertical  array  evenly 
spaced  in  depth  between  50  and  95  m.  Both  source  and  object 
are  midway  between  the  top  and  bottom  of  the  water  column. 
To  determine  the  presence  of  the  object  in  a  sound  channel,  we 
mode  form  as  the  object  passes  between  the  source  and  the 
receiver.  The  hydrophone  outputs  are  weighted  by  the 
amplitude  at  the  first  mode  and  summed.  This  array  response 
in  decibels  due  to  the  scattered  field  is  plotted  versus  distance 
along  the  track  in  Fig.  13.  The  maximum  response  usually 


comes  halfway  into  the  track  when  the  object  is  between  the 
source  and  the  receiving  array.  However,  since  the  scattered 
response  is  a  guided  wave,  it  can  take  a  dip  at  the  midpoint  as  it 
does  in  this  example.  The  array  response  due  to  the  scattered 
field  should  be  compared  with  the  array  response  due  to  the 
direct  field,  which  is  166  dB.  In  the  future  this  method  should 
prove  to  be  useful  in  discriminating  between  objects. 

V.  Conclusion 

We  have  examined  scattering  from  objects  in  an  unbounded 
environment,  in  an  environment  near  an  interface,  and  in  a 
waveguide.  Results  for  each  case  are  quite  different.  For  the 
unbounded  case  bistatic  and  monostatic  angular  distributions 
and  form  functions  are  important  tools  that  directly  relate  with 
objects.  For  objects  near  an  interface,  the  interaction  with  the 
incident  field  and  the  interface  is  a  very  important  factor.  For 
waveguide  scattering,  the  scattered  field  eventually  adapts  to 
the  waveguide  in  which  beamforming  techniques  would  be 
used  for  detection.  The  effect  of  motion  on  the  mode-formed 
or  beam-formed  array  then  becomes  crucial;  thus  transition 
time  becomes  important. 
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